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Segregation and mixing of granular mixtures during heap formation have important consequences 
in industry and agriculture. This research investigates three different final particle configurations 
of bi-disperse granular mixtures - stratified, segregated and mixed - during filling of quasi- two 
dimensional silos. We consider a larger number and relatively wider range of control parameters 
than previous studies, including particle size ratio, flow rate, system size and heap rise velocity. The 
boundary between stratified and unstratified states is primarily controlled by the two-dimensional 
flow rate, with the critical flow rate for the transition depending weakly on particle size ratio and 
flowing layer length. In contrast, the transition from segregated to mixed states is controlled by the 
rise velocity of the heap, a control parameter not previously considered. The critical rise velocity 
for the transition depends strongly on the particle size ratio. 



PACS numbers: 47.57.Gc, 81.05.Rm 

I. INTRODUCTION 

Heap flow of granular materials occurs in many con- 
texts [1, 2]. For example, when granular materials such as 
powders, grains, or pelletized polymers flow into the top 
of a container, a heap builds, where the granular material 
tumbles down the pile via a flowing layer that is a few 
particle diameters thick. When materials are mixtures of 
particles differing in size, density, shape, and/or surface 
properties, different components tend to distribute inho- 
mogeneously. In some cases, larger particles flow further 
down the heap than smaller particles resulting in segre- 
gation [3-10]. In other cases, large and small particles 
form alternating layers resulting in stratification [11-17]. 
In still other situations, large and small particles remain 
mixed [15]. 

In what were perhaps the earliest attempts to un- 
derstand the physical mechanisms driving segregation in 
heap flow, Williams [3, 4] and Drahun and Bridgwater [5] 
performed heap flow experiments using bi-disperse mix- 
tures of different-sized spherical particles. They proposed 
a percolation mechanism for heap segregation, in which 
small particles tend to sink through voids preferentially, 
while large particles rise to the free surface and roll to 
the end of the flowing layer. Consequently, small parti- 
cles accumulate below the upstream portion of the flow- 
ing layer at the center of the heap while large particles 
accumulate at the downstream end of the flowing layer of 
the heap adjacent to the bounding outer walls. Based on 
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FIG. 1. (a) Sketch of the top view of a 3D silo, where arrows 
indicate flow of material. Dashed box shows top view of the 
quasi-2D silo used here; (b) Side view sketch of a quasi-2D silo 
rising at the rise velocity v r = Q/A, where Q is volumetric 
feed rate and A = TW. The three stages of heap flow are: (I) 
Heap initiation; (II) Heap lateral growth; (III) Heap steady 
filling once it is constrained by the outer wall. 

this picture, Shinohara and co-workers [6, 7] developed a 
screening layer model based on conservation equations in- 
corporating the percolation mechanism. A recent exper- 
imental and computational study [10] performed to test 
the screening layer model [7] shows the model captures 
some key features of heap segregation in certain ranges 
of experimental parameters. However, some variables in 
the screening layer model (such as the penetration rate of 
segregating components and the velocity ratio of different 
sub-layers) can only be determined by fitting experimen- 
tal or simulation data, limiting the applicability of this 
model. 

Stratification in heap flow was first studied in detail 
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by Makse et al. [11, 12] 1 . In experiments, they found 
that when the components of granular mixtures differ in 
both size and shape, the mixtures stratify into alternat- 
ing layers of larger rough particles and smaller smooth 
particles. They attributed the formation of stratified lay- 
ers to competition between size segregation and shape 
segregation, and proposed that stratification only occurs 
for mixtures of large rough and small round particles. 
They also adapted a continuum model [12, 18, 19] and a 
cellular automaton model [11, 12, 19] to successfully re- 
produce the stratification observed in their experiments. 
However, Baxter et al. [15] found that stratification can 
also occur for different-sized smooth spherical particles, 
though they did not perform a systematic study over a 
wide range of parameters. Further, Baxter et al. [15] 
mentioned that a mixed state (no segregation or strati- 
fication) exists in certain situations such as at high feed 
flow rates. 



Although three final states of heap flow of bi-disperse 
granular materials (segregated, stratified, and mixed) 
have been observed and studied by different researchers, 
none of the past research appears to have investigated the 
dependence of the final particle distributions on a broad 
range of control parameters (see Table I), including vol- 
umetric feed rate Q, silo width W, 2D silo gap thickness 
T, species size ratio R = D\j D s (Di and D s are the large 
and small particle diameter, respectively), and absolute 
particle size. The focus of this research is to systemically 
explore how these parameters affect transitions between 
different final particle configurations. 

Figure 1(a) is a sketch of the top view of a three- 
dimensional (3D) silo, which typically comprises a ver- 
tical cylindrical container where granular material falls 
vertically along its centerline and flows radially in all di- 
rections down the heap, filling the container to its outer 
wall. Here, we use a quasi-two-dimensional (2D) silo, 
which can be thought of as a section of the 3D silo as 
shown in Fig. 1(a). This makes it easy to observe the 
final particle distributions and minimizes the volume of 
particles needed. In the quasi-2D silo, rather than feed 
rate Q, the relevant flow rate is the 2D volumetric flow 
rate down the slope at the peak of the heap, defined as 
q = Q/T, which decreases linearly along the flow direc- 
tion. As shown in Fig. 1(b), filling of a silo proceeds in 
three stages. In stage I, a somewhat irregularly-shaped 
initial heap forms. Shortly, the heap becomes angled in 
stage II and grows laterally until it reaches the bounding 
end wall. In stage III, the laterally constrained heap rises 
steadily at a constant rise velocity v r . We examine the 
steady filling stage (III), where the length of the flowing 
layer L is constant. 



One figure of Williams [4] shows stratification in heap flow, but 
the stratification was not noted. 



In this Article, we present an experimental study of 
heap segregation of granular mixtures of binary spheri- 
cal particles differing only in size in a quasi-2D silo for a 
larger number and relatively wider range of experimental 
parameters than in previous work (Table I). We observed 
all three final particle configurations - stratified, segre- 
gated and mixed. We find that for constant i?, the tran- 
sition from the stratified state to the unstratified state is 
controlled by q, while the transition from the segregated 
state to the mixed state is controlled by the rise velocity 
v r = q/W. These transitions depend on size ratio R. For 
the first time, phase diagrams are plotted to illustrate the 
effects of i?, W, and q or v r on transitions between differ- 
ent final particle configurations. We further provide new 
insight for the occurrence of the stratified configuration 
for granular mixtures of spherical particles compared to 
previous work [11, 15] and propose a dimensionless veloc- 
ity ratio - the ratio of rise velocity to particle percolation 
velocity - as a key control parameter for segregation. 

In the remainder of this paper, Sec. II describes the 
experimental setup. Sections III and IV present experi- 
mental results and discussion thereof. Section V presents 
our conclusions. 



II. EXPERIMENTAL SETUP 

The quasi-2D silo in our experiments consists of a pair 
of 91 cm x 69 cm x 1.27 cm vertical rectangular plates: 
one is thick plate glass for observation and measurement 
purposes and the other is aluminum to reduce electro- 
static charging. Vertical spacer bars were clamped be- 
tween the parallel plates to control the silo gap thickness 
T and to vary the silo width W . Granular mixtures were 
fed either at one end of the silo or at the centerline of 
the silo (the length of the flowing layer L in the latter 
case is half of the former case). Similar results were ob- 
served independent of whether the feed was at one end 
of the silo (total silo width W) or at the centerline of 
the silo (total silo width 2W). For the results presented 
here, values for W were 22 cm, 46 cm, 69 cm, and 91 
cm. We note that the size of the silo at W = 91 cm is 
comparable to small, full scale, industrial silos. The ef- 
fect of the silo gap thickness T on the transition between 
segregation and mixing is negligible, provided that T is 
more than four large particle diameters. However, strat- 
ification depends sensitively on T, as discussed in the 
Appendix. Here, we use T =1.27 cm for all experiments, 
unless otherwise noted. 

Seven combinations of different-sized soda-lime glass 
particles with size ratios R ranging from 1.3 to 6.0 were 
investigated (see Table II). We limit R to be less than 
(2/y/3 — 1) _1 = 6.464, above which spontaneous percola- 
tion occurs [20-22] (see Sec. IV C). To distinguish differ- 
ent species, different colors of particles were used. To en- 
sure roundness and similar surface properties such as fric- 
tion coefficient between different species, we purchased 
particles colored by the manufacturer (Sigmund Lind- 
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TABLE I. Experimental parameters in past and current research 
References Segregation type Size ratio Silo width 2D Silo thickness Flow rate Heap stage 







R 


W (cm) 


T (cm) 


Q (cm 3 /s) 


(see Fig. 1) 


Williams [3] 


size 


5.2 


31 a 


2.5 


unspecified 


II 


Drahun & Bridgwater [5] 


size 


1.3-2.0 


43 


unspecified 


unspecified 


III 


Shinohara et al. [6] 


size 


14.3 


30 


unspecified 


30-80 


II&III 


Shinohara et al. [7] 


size 


2.0-15.0 


15-18 


unspecified 


14-33 b 


III 


Baxter et al. [15] 


size 


2.0 


50 a 


9 


6, 736 c 


II 


Thomas [8] 


size 


2.0-50.0 


10 


3D silo 


27-40 a 


II 


Goyal & Tomassone [9] 


size 


1.3-5.0 


22 


0.5 


unspecified 


II 


Rahman et al. [10] 


size 


9.2-15.2 


15-18 


3D silo 


33 b 


III 


Makse et al. [11, 12] 


size & shape 


1.7-6.7 


30 


0.5-1 


unspecified 


II 


Grasselli & Herrmann [13] 


size & shape 


1.2-10.5 


30 


0.1-0.6 


0.2-3.5 c 


II 


Koeppe et al. [14] 


size & shape 


2.0 


27 


0.3-2.4 


0.5-7.4 c 


II 


Shimokawa & Ohta [16, 17] 


size & shape 


2.0-8.0 


60 


0.5 


0.08-1. 72 d 


II 


Current study 


size 


1.3-6.0 


22-91 


0.6-2.5 


1-420 


III 



a For particles fed at the silo center, W is half of the silo width. 

b Flow rate is estimated as Q = Q m /(0.6p m ), where Q m is mass flow rate and p m is material density. 
c Flow rate is estimated as Q = F/0.6, where F is feed flow rate on net particle volume basis. 
/(0.6pm), where Q m is mass flow rate and p m is estimated as 2.0 g/cm 3 . 



TABLE II. Binary mixtures of glass particles with equal mass 
fractions. 



TABLE III. Parameters 



Size ratio 
R 


Large particles 
Di (mm) a 


Small particles 
D s (mm) a 


1.3 


2.00±0.08 


1.51±0.09 


1.5 


1.69±0.05 


1.14±0.08 


2.0 


2.00±0.08 


1.00±0.06 


2.2 


1.10±0.06 


0.50=b0.05 


3.4 


1.69±0.05 


0.50±0.05 


4.0 


2.00±0.08 


0.50±0.05 


6.0 


2.98±0.05 


0.50=b0.05 



Parameter 



Mean particle diameter with standard deviation. 



ner GmbH, Germany). The granular mixtures in our ex- 
periments are composed of either metal-coated, surface- 
colored red and blue particles or clear and surface-colored 
black particles. Several trial experiments showed that fi- 
nal particle configurations are insensitive to the surface 
coatings of the two components. The material density of 
the particles is 2.59 g/cm 3 . 

A small auger feeder (Acrison Inc., NJ, USA) dispenses 
the granular mixtures into the silo. The auger feeder 
produces a stable and reproducible flow rate over a wide 
range of volumetric flow rates (1 to 420 cm 3 /s). Flow 
rate is varied by controlling the rotation frequency of 
the motor and the diameter of the auger. The granu- 
lar mixtures were composed of equal masses of large and 
small particles and were well-mixed upon filling the auger 
feeder. We performed several test experiments to mea- 
sure the mass fraction of the two species of different-sized 
particles after discharge from the auger feeder at several 
different flow rates and size ratios. We found that the 
mass fraction of the two components remained at 50:50 
(within ± 3%) during the entire experiment. Before each 



Description 



Type 



Q 
w 

T 

R — Di/D s 

q = Q/t 

v r q/W 



3D flow rate 

silo width 

silo gap thickness 

size ratio 

2D flow rate 

heap rise velocity 



controlled 
controlled 
controlled 
controlled 
calculated 
calculated 



experiment, anti-static spray (Sprayon, OH, USA) was 
applied to the glass wall to limit electrostatic effects. A 
digital camera in front of the glass wall recorded the fill- 
ing and the final state of the heap. We performed more 
than 400 experimental runs, systematically varying con- 
trol parameters (Table III) including the feed volumetric 
flow rate Q, system size (silo width W and silo thickness 
T), size ratio R, and absolute particle size. 



III. RESULTS 

Stratified, segregated and mixed final states noted by 
previous researchers were all observed in our experiments 
for stage III, as shown in Fig. 2. The thin free surface 
layer should be ignored as it is associated with resid- 
ual flow at the end of filling. At small flow rates, similar 
but more pronounced stratification than reported in Bax- 
ter et al. [15] occurs [see Fig. 2(a)]. The stratified state 
consists of alternating layers of small and large parti- 
cles parallel to the flow direction, coexisting with segre- 
gation along the flow direction, where the downstream 
region of the heap contains mostly large particles and 
the upstream region of the heap close to the feed point 
contains mostly small particles. Stratification becomes 



FIG. 2. (Color online) Three final particle configurations of different-sized granular mixtures in quasi-2D heap flow for R = 3.4, 
T — 1.27 cm, and W = 69 cm (the free surface layer is associated with residual flow at the end of filling and should be ignored): 
(a) Stratification at q — 0.8 cm 2 /s; (b) Segregation at q = 18.9 cm 2 /s; (c) Near mixing at q = 328 cm 2 /s. Dark (blue online): 
1.69 mm glass particles; light (red online): 0.5 mm glass particles. The width of each image is W. 



weaker and eventually disappears as q increases toward a 
critical value (to be discussed shortly). Above this crit- 
ical value full segregation is observed, as shown in Fig. 
2(b). In this regime, the heap consists of two distinct 
regions: The downstream region contains of nearly all 
large particles, while the upstream region consists of a 
few large particles scattered in a sea of small particles. 
The boundary between these two segregated regions is 
narrow. As q is further increased, the region of larger 
particles in the downstream portion of the heap shrinks 
and more large particles remain in the upstream portion 
of the heap. One might expect that at high enough q, 
the region containing all large particles disappears and a 
perfectly mixed state is achieved everywhere in the silo. 
However, due to limitations of our experimental appara- 
tus, the silo cannot always be filled at a high enough flow 
rate q to achieve perfect mixing. Instead, for most exper- 
iments at the highest achievable a near mixed state, 
see Fig. 2(c), is obtained where only a narrow large par- 
ticle region exists at the downstream end of the heap and 
the remainder of the heap is well-mixed. 



A. Stratification 

Stratification [see Fig. 2(a)] occurs at small flow rates 
over a wide range of size ratios and silo widths. We use 
the image intensity to quantify the final particle distribu- 
tions for stratification for each experimental run. Since 
the particles used in the experiments have different col- 
ors, the local image intensity, I, is monotonically related 
to local particle concentration. In all our experiments, 
small particles have higher intensity than large particles, 
so higher local intensity implies higher local concentra- 
tion of small particles. 

Figure 3 illustrates our method for measuring parti- 
cle concentration for stratification. We study the region 
outlined in Fig. 3(a), which is located in stage III of the 
heap formation and excludes the free surface region as- 
sociated with residual flow at the end of filling. The out- 
lined region is transformed to a rectangular box [see Fig. 
3(b)] by rotating by the angle of repose in the counter- 




15l - ■ — — ■ I 

90 100 110 120 130 

I (arbitrary units) 



FIG. 3. (Color online) Image processing method for quanti- 
fying stratification, (a) Image from the experiment with W = 
91 cm, T = 1.27 cm and q = 1.8 cm 2 /s showing a stratified 
mixture of 0.5 mm (light) and 1.1 mm (black) glass particles, 
(b) Dashed parallelogram region from (a) transformed to a 
rectangular box. (c) intensity I from (b) averaged over x and 
plotted as a function of y. 



clockwise direction and then "shearing" in the horizon- 
tal direction. The image intensity is averaged in the in- 
direction (0 < x < L) and plotted as a function of y [Fig. 
3(c)]. The intensity profile represents the variation of 
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FIG. 4. (Color online) ai/I showing decrease in stratification 
with increasing q for R — 2.2 at different W. Insets: images 
from experiments at the indicated data points, where D s — 
0.5 mm for small light particles, Di — 1.1 mm for large black 
particles, and W = 91 cm. 



species concentration due to stratification; the periodic 
intensity oscillation in the ^-direction corresponds to the 
alternating layers of small and large particles. 

To further quantify the stratification globally and to 
determine the transition to segregation at different ex- 
perimental conditions, we calculate the standard devia- 
tion of the intensity profile in Fig. 3(c) over the range of 

y coordinates, oj — yj^2^ =1 (li — I) 2 /(N — 1), where / is 
the mean intensity, Ii is the intensity at row z, and TV is 
the number of pixel rows in the ^/-direction. Larger ajjl 
indicates a higher degree of stratification. When there is 
no stratification, cji/I goes to a constant residual value 
of 0.005 associated with variations in lighting intensity 
and random fluctuations of the layer concentration. 

Figure 4 shows aj/I as a function of q for R = 2.2. 
c j / 1 is significantly larger at small q corresponding to 
strong stratification (long layers). As q increases, ai/I 
decreases as the stratified layers become shorter and 
stratification weakens. cfj/I decreases to a small con- 
stant value at a transitional 2D flow rate q t , where strat- 
ification disappears and only segregation occurs. The 
transition from a stratified state to an unstratified state 
occurs around q t = 6 cm 2 /s, independent of W. Similar 
trends are observed for all other R except the smallest 
value considered, R= 1.3, where no stratification occurs 
for all q and W tested. 

The influence of q , W, and R on the transition from 
stratified to unstratified states is examined by plotting 
a phase diagram as a function of q and i?, shown in 
Fig. 5. At each R, data for different values of W are 
artificially offset in three rows representing from top to 
bottom, W = 91 cm, W = 69 cm, and W = 46 cm, re- 
spectively, to show all of the data points. The three final 
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FIG. 5. (Color online) Phase diagram of final states (strat- 
ified: red O; segregated: black □; mixed: blue <0) in terms 
of q and R at three different W. Data are artificially offset 
in R to show each data point for different W: from top to 
bottom, W = 91 cm, 69 cm, and 46 cm, respectively. Hori- 
zontal solid lines denote the actual size ratio to guide the eye. 
Dashed lines segments mark the boundary between stratified 
and unstratified states. 



states are distinguished by symbols and the dashed line 
indicates the boundary between the stratified state and 
the unstratified state. (We discuss the transition from 
segregation to mixing in Sec. IIIB.) 

For 2 < R < 4, stratification occurs when q is less 
than 6 cm 2 /s. At these size ratios, the transition from 
the stratified state to the unstratified state occurs at the 
same transitional flow rate q t for all three values of W. 
For R = 6, stratification is observed up to q = 12 cm 2 /s 
for all three values of W. This increase in q t is possibly 
because R is close to the size ratio for spontaneous seg- 
regation, 6.464, so that other factors such as wall effects 
or horizontal segregation due to spanwise shear rate gra- 
dients [23] may significantly influence the stratification. 

At small i?, stratification diminishes. When R =1.3, 
no stratification occurs for all W. When R = 1.5, strati- 
fication occurs only at W = 69 cm and 91 cm; no strat- 
ification is observed at W = 46 cm. Thus, W may also 
affect the occurrence of the stratification at small R as 
discussed further in Sec. IV A. 



B. Segregation 

Full segregation occurs when q increases beyond the 
transitional value q t . A similar image processing method 
to that for stratification is used to quantify segregation. 
As shown in Fig. 6(a), the outlined region in stage III is 
considered. The parallelogram is transformed to a rect- 
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(a) 200 




FIG. 6. (Color online) Image processing method for quan- 
tifying segregation, (a) An image from the experiment with 
W = 69 cm, T — 1.27 cm, and q — 10.9 cm 2 /s for a mixture 
of 0.5 mm (light) and 1.1 mm (black) glass particles, (b) The 
region in the dashed parallelogram from (a) transformed into 
a rectangle, (c) intensity / from (b) averaged over y and plot- 
ted as a function of x/L. AL denotes the width of the dark 
region of large particles at the end of the heap, and L is the 
length of the flowing layer. 



angle as in Fig. 6(b), and the image intensity is averaged 
in the ^-direction and plotted as a function of x as shown 
in Fig. 6(c) for a typical case. The upstream portion of 
the heap (x/L < 0.6) has a smoothly varying concentra- 
tion of particles of the two different sizes, whereas the 
downstream end (x/L > 0.6) has a nearly uniform con- 
centration of only large particles. 

Figure 7(a) shows a series of intensity profiles plotted 
as a function of x/L at different W and q at R = 2.2. The 
profiles overlay one another for different values of W, as 
discussed shortly. Each profile represents the final state 
distribution of the two segregated species. The concen- 
tration of small particles is higher in the upstream region 



W = 46 cm 
W = 69 cm 
W = 91 cm 




200 
100 




vr = 1.2 cm/s 

) 6.2 ~~ OA 06~~ ~~ 08 

x/L 




• W = 46 cm 

♦ W = 69 cm 
■ W = 91 cm 




1 2 3 4 5 

v r (cm/s) 

FIG. 7. (Color online) (a) Intensity I as a function of x/L 
for a mixture of 0.5 mm (light) and 1.1 mm (black) particles 
at v r = 0.08 cm/s (top), 0.24 cm/s (middle) and 1.2 cm/s 
(bottom), respectively, (b) AL/L vs. v r for the same mixture 
at different silo widths W. Insets show images for W = 69 
cm and indicated data points. Dashed line is the fit AL/L — 
(v r /a)~ h , where a =0.017 cm/s and b =0.493. 



and the concentration of large particles is higher in the 
downstream region. The boundary between these two re- 
gions is narrow (less than 0.2L). Close to the feed zone 
of the heap (x/L < 0.2), the concentration of small par- 
ticles is slightly smaller than in the rest of the upstream 
region [Such as Fig. 6(c) and Fig. 7(a)]. This likely occur 
because the incoming particles discharged from the auger 
feeder start flowing from a nearly stationary state after 
falling onto the heap so that the segregation is weaker in 
this region. 

As indicated in Fig. 7(a), the similar intensity distribu- 
tions in each plot correspond to identical rise velocities, 
v r , but different silo widths W and flow rate q. In stage 
III, v r = Q/(WT) = q/W and is independent of the an- 
gle of repose that changes slightly as q or R varies [24]. 
Thus, when W and q are varied together, v r , instead 
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FIG. 8. (Color online) Phase diagram of final state (strat- 
ified: red O; segregated: black □; mixed: blue 0) in terms 
of rise velocity v r and size ratio R at three silo widths W. 
Data are artificially offset in R to show each data point for 
different W: from top to bottom, W — 91 cm, 69 cm, and 
46 cm, respectively. Horizontal solid lines denote the actual 
size ratio to guide the eye. Dashed line segments indicate the 
boundary between segregation and mixing. 

of controls the final particle distributions for segrega- 
tion. In other words, for a given mixture at the same v r , 
small and large particles distribute similarly along the 
flow direction for different values of W and q. When v r 
increases, more large particles stay in the upstream re- 
gion so that the concentration of small particles decreases 
and the width of the downstream region of large particles 
AL decreases as well, as shown in Fig. 7(a). However, the 
concentration of large particles in the downstream region 
of the heap does not change as v r increases (the intensity 
at large x/L is constant at different v r ). Based on our 
experiments, there is almost always a region at the end 
of the heap close to the end walls containing nearly all 
large particles, regardless of the flow rate, silo width, and 
size ratio 2 . 

To determine the transition from segregation to mix- 
ing as a function of v r and R, we quantify the degree 
of segregation using the normalized length of the large 
particle region AL/L. This measure is similar to the ap- 
proach of Shinohara et al. [6], who used 1 — AL/L. In 
stage III, AL/L represents the approximate mass of seg- 



2 There are two exceptions: (1) when R < 1.5, the large particle 
region at the end of the heap is absent in many cases, which 
is discussed in Sec. IV C. (2) At large D s or small W, a narrow 
region of small particles appears at the end of heap due to another 
mechanism: the bouncing of particles. For all data presented in 
this paper, the bouncing-induced segregation is minimized by 
using small particles and large W (see Sec. IV D). 



regated large particles relative to the mass of the mixture. 
Therefore, based on mass conservation, < AL/L < 0.5, 
where AL/L — means perfect mixing and AL/L = 0.5 
means complete segregation. 

Figure 7(b) plots AL/L as a function of v r for different 
W at R = 2.2. At other size ratios similar trends occur. 
As v r is increased, AL/L decreases, indicating that seg- 
regation becomes weaker at higher v r . For v r > 1.5 cm/s, 
AL/L < 0.1, which means more than 80% of the large 
particles are mixed with small particles in the upstream 
region of the heap (x/L < 0.9). This corresponds to a 
mass ratio of small particles in this region of 56%. Curves 
of AL/L as a function of v r at different W collapse onto 
a single curve, consistent with the collapse shown in Fig. 
7(a) for different W. 

Due to limitations of our experimental apparatus (par- 
ticles overflowing the entrance of the silo), a rise velocity 
greater than 10 cm/s (equivalent to filling the entire silo 
in less than 6 s) cannot be achieved. A perfect mixed 
state is therefor difficult to obtain at most R except for 
R = 1.3, where relatively good mixing is observed for 
v r > 1 cm/s. We fit the experimental data in Fig. 7(b) to 
a power law as AL/L = (v r /a)~ b using the least squares 
method for each size ratio [see Fig. 7(b) for R = 2.2], 
where the fitting parameters a and b depend on R. Based 
on this fit, a cutoff value for AL/L = 0.15 is selected to 
determine the transitional rise velocity between segrega- 
tion and mixing, so that the boundary can be determined 
at different R. This value corresponds to a small particle 
mass ratio of 59% in the region x/L < 0.85. 

Using AL/L = 0.15 as the cutoff value between seg- 
regation and mixing, a phase diagram similar to that in 
Fig. 5, but in terms of v r instead of q, is constructed 
to quantify the effects of i?, W and v r on the transition 
from segregation to mixing in Fig. 8. The dashed line 
segments in the figure marks the approximate bound- 
ary between segregation and mixing. The phase diagram 
demonstrates that segregation transitions to mixing at 
the same v r for different W at the same R. v r at the 
transition increases by roughly one order of magnitude 
as R is increased from 1.3 to 6.0. 



IV. DISCUSSION 

As shown in Sec. Ill, stratification generally occurs at 
low flow rates, and the transition from stratified to un- 
stratified states depends on the flow rate q. In contrast, 
the transition from segregated to mixed states depends 
on the rise velocity of the heap v r . In this section, we dis- 
cuss possible mechanisms for these phenomena. Further, 
we briefly discuss, with respect to segregation, the lim- 
its of size ratio and the effect of particle bouncing after 
initial impact with the heap. 
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FIG. 9. (Color online) Static angle of repose for monodis- 
perse glass particles at different silo gap thicknesses. Results 
are for 0.5 mm, 1.1 mm, 1.7 mm, 2 mm, and 3 mm monodis- 
perse particles flowing into the silo at identical feed rate. The 
static angle of repose is measured after filling is stopped. 

A. Stratification dynamics 

As mentioned in Sec. I, stratification has been most of- 
ten observed in mixtures of large rough and small smooth 
particles as reported by Makse et al. [11, 12]. The strat- 
ification is attributed to differences in the angle of re- 
pose: steeper for rough particles than for smooth spheri- 
cal particles. Although Williams [3] and Baxter et al. [15] 
showed some evidence of stratification for different-sized 
smooth spherical particles, the existence of stratification 
in mixtures of different-sized spherical particles has been 
debated. 3 Furthermore, these works [3, 15] only provide 
a few limited examples of stratification. 

In this research, particles are all spherical so shape 
effects are excluded. Furthermore, except for the smallest 
particles, the angle of repose is nearly independent of 
particle size at the same silo gap thickness T, as shown 
in Fig. 9. As described in Sec. Ill A, stratification occurs 
for bi-disperse spherical particles over a wide range of 
flow rates, size ratios, and system sizes, which indicates a 
different mechanism for stratification than that proposed 
by Makse et al. [11, 12]. 

Careful observation of stratification for spherical parti- 
cles in our experiments indicates that the appearance of 
layers of particles of different sizes is associated with the 
occurrence of discrete avalanches at low flow rates. As 
Fig. 10 shows, during an avalanche, large (dark) parti- 
cles segregate via percolation to the free surface and roll 



3 Though Williams [3] and Baxter et al. [15] observed stratification 
with smooth spherical particles, Makse et al. [11, 19] argued that 
some shape induced segregation may have occurred. 




FIG. 10. (Color online) Illustration of stratification forma- 
tion at R = 3.4, W = 91 cm and q — 0.9 cm 2 /s. (a) and 
(b) show consecutive occurrences of two avalanches. Kink 
A is formed in the first avalanche, and it moves downstream 
when the second avalanche occurs, and during which kink B is 
formed. Dark (blue online): large particles; light (red online): 
small particles. 



down the free surface forming the front of an avalanche. 
As time passes, the large particles in the front stop flow- 
ing and form a subtle kink in the slope of the heap [kink 
A in Fig. 10(a)]. Upstream from the kink, a stripe made 
of a thick layer of small particles at the bottom and a 
mono-layer of large particles at the surface forms [Fig. 
10(a)]. Since the local angle of repose at the kink is larger 
than the static angle of repose, the kink is metastable. 
When a new avalanche occurs upstream of the kink and 
on top of the previous avalanche, kink (A) propagates 
further downstream, resulting in a part of the stripe as- 
sociated with kink A reflowing with the arrival of the new 
avalanche [kink B in Fig. 10(b)]. As a result, the layer 
of small particles in the original stripe moves downward 
until the front of the avalanche, which consists of large 
particles, reaches the end wall of the silo, forming a layer 
of small particles above the larger particles from the pre- 
vious kink (A) and below the larger particles from the 
current kink (B). In contrast, at high flow rates, surface 
flow is continuous, and large particles do not form kinks 
since there are no avalanches. Instead, large particles 
are continually advected to the downstream region of the 
heap. Small particles settle in the upstream region of the 
heap, resulting in the typical segregation pattern shown 
in Fig. 2(b). 



In addition, stratification is also limited by size ratio, 
silo width, and silo gap thickness. At small size ratios 
(e.g. R = 1.3), segregation in the flowing layer evolves 
slowly that the initial stripe of small and large particles 
is weak. The silo also needs to be wide enough so that a 
kink can form somewhere along the slope of the heap. For 
example, no stratification occurs for R = 1.5 and W = 46 
cm, while stratification does occur for large widths W 
(see Fig. 5). Silo gap thickness also affects the degree 
of stratification as discussed in the Appendix. A more 
detailed study of stratification of spherical particles is 
underway. 



B. Control parameters for segregation 

The dynamics and mechanisms for segregation in heap 
flow have been previously studied [3-7] as discussed in 
Sec. I. At its simplest, the competition between percola- 
tion of two different-sized components perpendicular to 
the flow direction and the advection of the mean flow in 
the flow direction determines the degree of final segrega- 
tion. For instance, if percolation of small particles down- 
ward through the flowing layer takes longer than the time 
to than the time to reach the downstream of the heap, 
more small particles will accumulate at the downstream 
region of the heap along with large particles resulting in a 
more mixed state. On the other hand, when small parti- 
cles percolate quickly to the bottom of the flowing layer, 
they remain in the upstream region of the heap, resulting 
in stronger segregation. 

Based on this picture, we define a dimensionless time 
(t) = tp/td, where t p represents the time scale for perco- 
lation, and td represents the time scale for downstream 
convection by the mean flow. Based on this definition, 
(t) ^> 1 indicates better mixing, while (t) <C 1 indi- 
cates stronger segregation. Assuming a constant average 
downstream velocity Vd and linear velocity profiles in the 
flowing layer for simplicity (similar to previous work on 
mono-disperse heap flow [18, 25]), we obtain td = L/vd, 
t p = S/vp and q = UdS = v r L (based on mass conserva- 
tion in the flowing layer), where 5 is the mean thickness 
of the flowing layer and v p is the mean percolation ve- 
locity. With these assumptions, the dimensionless time 
scale can be expressed as a velocity ratio 



(t) V r /Vp 



(1) 



Equation (1) shows that for a particular percolation 
velocity v pi which depends on the size ratio i?, the degree 
of final segregation depends only on v r . Increasing v r 
increases (£), indicating a decrease in segregation and 
an increase in mixing. For constant v r , (t) decreases as 
v p increases (e.g. R increases). This means that at the 
same flow conditions, mixtures with large size differences 
should have a higher degree of segregation. Figure 11(a) 
shows profiles of AL / L as a function of v r for several size 
ratios (R = 1.5,2.2,3.4,4,6). As predicted by Eq. (1), 




0.5 1 

<t> ~ v r / Vg(R-1)D s 

FIG. 11. (Color online) (a) AL/L vs. v r at different R from 
1.5 to 6. At each R, data for different values of W use the 
same symbols, (b) AL/L as a function of the dimensionless 
time scale (£) ~ v r j \J g{R — 1)D S . 



AL/L, a measure of segregation, decreases as v r increases 
at constant R. Furthermore, at constant v r , as v p is 
increased by increasing R, AL/L increases, indicating 
stronger segregation at larger R. 

To further examine the dependence of segregation on 
(t) while varying v r and v p simultaneously, a relation be- 
tween the percolation velocity v p and experimental pa- 
rameters such as R and strain rate at different q is needed. 
Bridgwater and colleagues [26, 27] systemically studied 
the influence of different parameters including the size 
ratio, density ratio, strain rate, and normal stress, on 
percolation velocities in various sheared systems. They 
found that particle size ratio has the greatest influence 
on the percolation velocity, while the others have little 
effect. However, an analytic relation is lacking. Here, 
to first order of approximation, we assume v p depends 
only on R along with gravity, which drives the flow. 
Figure 11(b) shows AL/L as a function of (t) assuming 
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FIG. 12. (Color online) Intensity as a function of x/L for 
R = 1.3 at q = 1.8 cm 2 /s (black 0) and g = 21.6 cm 2 /s (red 
O) for (a) W = 46 cm and (b) W = 91 cm. 

v p oc y/g(R — 1)D 3 , where — 1 is used to assure zero 
percolation velocity for monodisperse particles. Curves 
of AL/L at different R collapse well compared with Fig. 
11(a). However, for R = 1.5, the scaling does not col- 
lapse the curves for (t) > 0.07 (corresponding to v r > 0.5 
cm/s), presumably because in this situation, other mech- 
anisms such as ordinary diffusion also play important 
roles in segregation. Thus, the expression we use for the 
percolation velocity is not always applicable. 



C. Upper and lower limits of size ratio 

In this work we considered size ratios from 1.3 to 6.0, a 
relatively large range compared to previous studies (Ta- 
ble I) of bidiperse mixtures of spherical particles. We 
choose R < 6 for two reasons. First, as shown in Refs. 
[20-22] and references therein, when R is larger than 
(2/>/3 - l)" 1 = 6.464, a small particle can percolate 
through the smallest voids between three large particles 
without external agitation. Thus, the mechanism for seg- 
regation is different from that at the smaller R investi- 
gated here. Further, as mentioned by Fan and Hill [23], 
a shear gradient can also drive segregation in the span- 
wise direction across the silo. At larger R , the shear 
effect is more important than at smaller R. During the 
experiments, at R = 6, we observed signs of both effects. 
Small particles preferentially fill the voids between large 
particles at the side walls, which causes a variation of 
particle concentration between the side walls. Based on 
observations at the top of the free surface, large particles 
also tended to gather toward the middle of the gap away 
from the side walls for R = 6, probably due to the shear 
effect [23]. 

At small size ratios (e.g. R = 1.3), no stratification is 



observed because of the small size difference, as discussed 
in Sec. IV A. Segregation also becomes much weaker at 
small size ratios and exhibits different characteristics 
than at other size ratios. As mentioned by Goyal and 
Tomassone [9], when R is smaller than a critical value 
of 1.4, the concentration of each species changes contin- 
uously along the heap, as shown for small q in Fig. 12(a) 
for R = 1.3. For R > 1.4, segregation is more complete 
and boundaries between the two segregated regions are 
sharper, as shown in Fig. 7. Our results show that the 
transition between these two different segregated states 
additionally depends on q and W, as shown in Fig. 12. 
At W = 46 cm for R = 1.3 and within the entire range of 
q in our experiments, continuously varying segregation is 
always observed [Fig. 12(a)]. However, for W = 91 cm, 
more complete segregation occurs at lower q and continu- 
ously varying segregation occurs at higher q [Fig. 12(b)]. 
One possible reason for this behavior is that at small 
i?, ordinary diffusion, which causes re-mixing of different 
components, may become comparable in importance to 
percolation over shorter W. Diffusion in flowing gran- 
ular materials is related to the flow kinematics such as 
local shear rates [28], which are closely associated with 
the flow rate and system size. At small flow rates or 
large silo widths, diffusion could be weaker than perco- 
lation so that segregation is stronger. Further study of 
the segregation mechanism at small size ratios R < 1.5 
is currently underway. 



D. Bouncing-induced segregation 

When the falling particles feeding the heap impact the 
top of the heap, they sometimes bounce away from the 
top of the heap toward the end of the silo. Since small 
particles generally gain momentum when colliding with 
large particles (similar to when a lighter golf ball col- 
lides with a heavier basketball), they tend to bounce fur- 
ther down the heap than large particles. This can result 
in segregation that is reversed and unrelated to that in 
the flowing layer. Bouncing-induced segregation, briefly 
mentioned by Drahun and Bridgwater [5], causes small 
particles to segregate to the downstream region of the 
heap, which is, of course, opposite to what occurs during 
the free surface segregation studied in this paper. 

Several factors can influence bouncing-induced segre- 
gation including the silo width, the fall height of the par- 
ticle feed stream, and the relative and absolute particle 
size. Figure 13 shows the influence of these factors on 
bouncing-induced segregation. In Figs. 13 (a)- (c), W is 
varied while all other parameters are fixed. For the small- 
est silo width, W = 22 cm [Fig. 13(a)], a narrow vertical 
band of small light-colored particles forms at the down- 
stream end of the heap (right wall). This is opposite to 
what is observed under normal free surface segregation 
conditions. When W is increased to 46 cm and 69 cm 
[Figs. 13(b)-(c)], no small particles are observed adjacent 
to the right end wall. Even though bouncing of particles 
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FIG. 13. (a)-(c) Images showing reduction in bouncing- 
induced segregation with increasing silo width W at R = 4 
with 2 mm (black) and 0.5 mm (light) glass particles and 
q = 14 cm 2 /s. (a) W = 22 cm, (b) W = 46 cm, and (c) 
W = 69 cm. (d) and (e) Images at similar size ratio (R « 
2) showing reduction in bouncing-induced segregation with 
decreasing absolute particle size at W = 69 cm and q = 140 
cm 2 /s. (d) 2.0 mm (black) and 1.0 mm (light) particles, (e) 
1.1 mm (black) and 0.5 mm (light) particles. 



still occurs at larger W, the bouncing particles re-enter 
the flowing layer before reaching the end of the silo, so 
that bouncing-induced segregation is negligible. The ef- 
fect of fall height on bounce-induced segregation can be 
seen by considering a single experimental run [see Fig. 
13(a)]. As the heap grows, the fall height decreases, so 
that the bouncing-induced segregation becomes weaker. 
Consequently, in Fig. 13(a), the small particle region at 
the right wall of the silo becomes narrower from bottom 
to top. At the smallest fall height near the end of the 
experiment and corresponding to the top of the vertical 
band of particles, there are only large particles close to 
the end wall, indicating a strong decrease in bouncing- 
induced segregation. In Figs. 13(b)-(c), in addition to the 
fact that no small particles are observed in the down- 
stream regions of the silo next to the right wall, the 
boundaries between large and small particle region are 
vertical (essentially independent of the fall height), indi- 
cating that bouncing-induced segregation has negligible 
influence on the final particle distributions. 

Figures 13(d)-(e) compare segregation when the granu- 
lar mixtures have similar size ratios (R « 2) but different 
absolute sizes [2 mm and 1 mm in Fig. 13(d) and 1.1 mm 
and 0.5 mm in Fig. 13(e)]. The mixture with the larger 
absolute size exhibits stronger bouncing-induced segrega- 
tion [small particles accumulate at the downstream wall 
of the silo in Fig. 13(d)] than the mixture of smaller ab- 
solute size [no small particles accumulate at the down- 
stream wall of the silo in Fig. 13(e)]. The inertial and 
gravitational forces of the particles are proportional to 
the cube of particle radius, while the air drag is at most 
proportional to the square of particle radius (at high 



Reynolds numbers) [29]. Therefore, the ratio of inertial 
or gravitational forces to air resistance is greater for the 
mixture of larger particles than for the mixture of smaller 
particles. As a result, the small particles in the mixture 
bounce farther after impact [Fig. 13(d)]. 

This cursory study of bouncing-induced segregation 
demonstrates that bouncing-induced segregation can 
compete with free surface segregation and may result in 
different final segregation patterns. However, as long as 
the silo width is large enough, bouncing-induced segre- 
gation is minimized and free surface segregation domi- 
nates the final particle distributions. To further study 
the effects of bouncing-induced segregation in heap flow, 
alternative methods to feed particles to the top of the 
heap (e.g. similar to that in Ref. [30]) may be needed 
to isolate the effects of these two different segregation 
mechanisms.. 



V. CONCLUSIONS 

In this paper, we have shown that three different final 
configurations in heap flow - stratified, segregated and 
mixed - can be obtained by controlling flow properties, 
either the flow rate q or the rise velocity of the heap 
v r . Stratification is associated with discrete avalanches 
at low flow rates (q smaller than ~10 cm 2 /s for spherical 
glass particle mixtures) at most size ratios (R > 1.4). 
The silo width W or alternatively the flowing layer 
length L, should also be large enough that a kink can 
form along the slope of the heap (instead of at the end 
of the heap) during the avalanche. The transition from a 
stratified to an unstratified state is governed by q. When 
q is larger than a transitional value dependent on i?, the 
heap flow transitions from a discrete avalanche regime 
to a continuous flow regime with segregation in which 
neither kinks nor stratification are observed. The degree 
of segregation is determined by competition between ad- 
vection by the mean flow and percolation through the 
flowing layer. This competition can be characterized by 
the ratio of the rise velocity of the heap v r to the perco- 
lation velocity v p , which mainly depends on R. At the 
same i?, as v r increases, the degree of segregation de- 
creases and eventually transitions to a mixed state. The 
transitional rise velocity becomes larger as R increases. 

There are three major points, among others, that this 
study raises: (i) Stratification of different-sized spherical 
particles is observed for a wide range of flow rates and 
size ratios, but the dynamics of stratification appear dif- 
ferent from those for stratification of different size and 
shape particles observed in Makse et al. [11, 12]. These 
apparent differences in the physical mechanisms for strat- 
ification need further study; (ii) This study shows that 
other mechanisms play important roles in heap segrega- 
tion for the following situations: a) At small size ratios 
(R <1.5), ordinary diffusion appears to become compa- 
rable to percolation, suggesting that a different segre- 
gation configuration occurs (continuously varying segre- 
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FIG. 14. (Color online) AL/L vs. v r for mixtures of 2.98 mm 
and 0.5 mm particles (R = 6) at different silo gap thicknesses 
T and silo widths W, showing segregation is insensitive to T 
forT>4A- 



gation). b) For R near the large size ratios at which 
spontaneous percolation can occur (R « 6.5), wall ef- 
fects or shear-induced segregation may cause spanwise 
segregation between the side walls, which requires study 
of segregation/stratification in the bulk of a fully 3D silo, 
(iii) When particles have sufficient impact velocity or the 
silo is not wide enough, bouncing of particles after heap 
impact can cause segregation opposite that of the usual 
segregation in the flowing layer. 

We have only studied quasi-2D configurations of heap 
flow. While we expect that these results are indicative of 
those in 3D heaps, further work is needed to confirm this. 
Furthermore, we have only considered the steady filling 
stage of heap formation [see Fig. 1(b)], leaving initial 



heap formation and heap growth for future investigation. 
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Appendix: Effects of quasi-2D silo gap thickness T 

The effect of quasi-2D silo gap thickness T on stratifi- 
cation and segregation was considered for several binary 
mixtures at various q to assure that T has no significant 
influence on the results. Silo thickness T was varied from 
0.64 cm to 2.54 cm (corresponding to 3 < T/D t < 12) 
and q was varied from 1.2 cm 2 /s to 331 cm 2 /s. Measure- 
ments of stratification and segregation are carried out 
using the same criteria as those in Sec. III. 

The effects of T on segregation are shown in Fig. 14 
by plotting AL/L vs. v r for R = 6. Within the range 
of T investigated, the degree of segregation does not de- 
pend on T, as long as T > 4Di. The same trends are 
observed for other values of R. In contrast, stratification 
is influenced by T, but no clear trend is observed. For 
example at R = 6, when T increases from 1.27 cm to 
2.54 cm, a /I decreases at the same q. At T — 2.54 cm, 
a/I is fairly small over the entire range of indicating 
that at this T, stratification is barely observable for all 
experimental parameters. However, at R = 2 or R = 3.4, 
and certain values of a /I at larger T is larger than at 
smaller T, indicating that stratification does not always 
monotonically decrease when T increases at all R. The 
mechanism for the dependence of stratification on T re- 
mains unclear and needs further investigation. In this 
paper, we use T = 1.27 cm for all experiments as it is 
sufficient for achieving T-independent segregation and it 
minimizes the volume of particles needed to perform the 
experiments. 
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